Fast transport of Bose-Einstein condensates 
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We propose an inverse method to accelerate without final excitation the adiabatic transport of 
a Bose Einstein condensate. The method, applicable to arbitrary potential traps, is based on a 
partial extension of the Lewis- Riesenf eld invariants, and provides transport protocols that satisfy 
exactly the no-excitation conditions without constraints or approximations. This inverse method is 
complemented by optimizing the trap trajectory with respect to different physical criteria and by 
studying the effect of noise. 
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Introduction. — A major goal of atomic physics is the 
comprehensive control of the atomic quantum state for 
fundamental research and applications in interferome- 
try, metrology, or information processing. The ability 
to manipulate Bose-Einstein condensates may be par- 
ticularly rewarding due, for example, to their potential 
in interferometric sensors, but it is also challenging, as 
their low temperatures make them more fragile than or- 
dinary cold atoms. A basic operation is the transport of 
the condensate to appropriate locations such as a "sci- 
ence chamber" , or to launch or stop the atomic cloud. 
This transport has been performed with several tech- 
niques based on adiabatic, slow motion to avoid exci- 
tations and losses [1,-^]. Long transport times, however, 
may be counterproductive since the condensate is more 
exposed to noise and decoherence, and also limit severely 
the repetition rates and signal to noise ratios. Fast, non- 
adiabatic but "faithful" transport of cold atoms, i.e., 
leading to the desired final state, has also been inves- 
tigated experimentally [5] and theoretically For 
the Schrodinger equation (SE) an inverse engineering 
method based on constructing Lewis-Riesenfeld invari- 
ants and corresponding dynamical modes (solutions of 
the SE formed by invariant eigenvectors times a phase 
factor) provides a "shortcut to adiabaticity" How- 
ever the invariant concept, i.e., an operator satisfying 
dl/dt = dl{t)/dt+ [I{t),H{t)]/{ih) = with constant 
expectation values for arbitrary states that evolve with 
the Hamiltonian H , is not directly applicable to the non- 
linear Gross-Pitaevskii equation (GPE). In fact previous 
extensions of this inverse technique to expansions of con- 
densates required special regim es or time-dependent Fes- 



hbach resonance control [llHlSj. By contrast, we show in 



this letter that the transport of a condensate scales in the 
right way for applying a generalized inverse method, so 
that no approximations or limits are required to design 
fast processes without final excitation. Strictly speaking 
we do not construct invariants as in the linear theory 



but, as far as ground-state to ground-state transport is 
concerned, the formal structure of the dynamical trans- 
port modes of the linear case remains valid for the GPE. 
This will be illustrated with a numerical example and 
compared with a direct method. In the final part of the 
article we shall optimize the trap trajectory according 
to several criteria, since the invariant-based inverse engi- 
neering provides a family of possible transport solutions, 
an analyze the effect of noise. 

General setting. — Our starting point is the GPE for po- 
tentials whose Schrodinger dynamics admit a quadratic 
invariant in momentum [8, 14, isj ]. 
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where f = f{t) is arbitrary, Vq is the Laplacian in Garte- 
sian coordinates for £) = 1, 2, or 3 dimensions, U is an ar- 
bitrary potential function of the argument cr = {q—a)/p 
and a;(t), the force F{t), a — a{t) and the scaling func- 
tion p = p{t) satisfy 
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p + u\t)p, 
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where wq is a constant and the dots represent time deriva- 
tives. The physical meaning of a. depends on the pro- 
cess type, see [8|, and will be clarified below. The wave 
function is normalized to one. Without the non-linear 
term, an arbitrary solution of Eq. ([T]) can be written as 
a linear combination of the eigenvectors ipn of the dy- 
namical invariant / [31, il){c{,t) — c„e*""^*^?A„(q, t), 
I{t)tpn{{l,t) = A„V'n(q, i)j where the amplitudes c„ and 
the eigenvalues An are constants. Here ipn is normal- 
ized to one, but continuum, delta-normalized states are 
also possible. The Lewis-Riesenfeld phases an satisfy 
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^ ^ = ("^n \^^m ~ ^\ V'") 0' A single mode so- 
lution of the SE = 0) takes the form 
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where we have introduced a scaled time r = J^dt'p~^, 
and (/)„ satisfies a Schrodinger equation with a time- 
independent Hamiltonian. 

These results can be generalized partially for the GPE, 
and the extent of the generalization depends on the pro- 
cess type. Inserting Eq. (|4]) as an ansatz for a time 
dependent solution of the GPE, 0„ must satisfy 
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Unlike the linear theory, we cannot construct the general 
solution by linear superposition, so we restrict the treat- 
ment to a single mode, for example the ground state, 
and therefore we drop the n subindex in Eq. ([5]) and 
hereafter. 

Equation ([5]) is very general and applicable to com- 
pressions, expansions, or transport for harmonic or an- 
harmonic potentials. It is most useful when p^^^g^ does 
not depend on time, since the physical solution of the 
time-dependent problem is then mapped, via Eq. (j4]) , to 
the solution of a much simpler stationary equation. This 
happens in several physically relevant cases, in particular 
for expansions of BECs when D = 2, or by tuning g^^ as a 
time-dependent coupling to cancel the time dependence 
of p^~^ Different time scalings combined with a 

Thomas-Fermi approximation also lead to a stationary 
equation 

Transport Processes. — Here we are interested in the 
simple but very important case p(t) = 1 Vt, associated 
with rigid transport processes. Then t = t and the coef- 
ficients of Eq. ([5]) are time independent. It is also useful 
to define (f>{(T,t) — e~*^*/''x(''')i where p is the chemical 
potential and xl"") satisfies the stationary GPE 
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The physical solution of the time-dependent GPE equa- 
tion for a single transport mode is 

(7) 

which is a fundamental result. Since this wave function is 
shape invariant the only possible excitations associated 
with such a mode are center of mass oscillations with 



constant mean field energy. In the following we shall 
apply it to ID transport and omit vector notation. 

Inverse engineering by harmonic transport. — In a ID, 
horizontal, harmonic transport of a condensate from to 
d in a time tf with zero mean velocity at i = and tj, 
a (now a scalar function) must be chosen to match the 
transport mode ([7]) with the instantaneous eigenstates of 
the Hamiltonian, including the mean field term, at times 
t — and t — tf. 

A concrete example will illustrate how this works. 
Consider a ^^Rb BEG in the F = 2, = 2 ground 
state constituted by 3000 atoms The transport of 
BECs aided by microchips can use a "bucket chain" ISl , 
or a single harmonic and frequency-stable bucket 19!]. 
We assume here a single bucket with = 27r • 50 Hz 
moved from go(0) = at time i = to qo{tf) = d = 1.6 
mm at tf. The time-dependent GPE is now 
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a particular case of Eq. ([T]) with oj{t) = wq, C/ = 0, 
F{t) = mw^qoW, f{t) = mujlql{t)/2 and p(t) = 1, so 
Eq. ([2]) does not play any role and a = gc has to satisfy 



qc + i^l{qc - go) = 0, 



(9) 



the equation for a classical trajectory qc{t) in a moving 
harmonic potential. (Note that for an abrupt shift of the 
trap one recovers the scaling for dipole oscillations |2^.) 
If we impose at < = the initial conditions 



qM = 9c (0) = qM = 0, 



(10) 



the transport mode ([7]) becomes equal to the instanta- 
neous eigenstates of Eq. ([8]) at i = 0. To solve Eq. ([8]) 
we proceed in two different ways, using direct or inverse 
approaches. 

In the direct approach we fix first the evolution of the 
center of the trap qo(i)- for example, see Fig. 5 

there, qo(t) is increased linearly during a quarter of the 
transported distance (i/4, then kept constant for (i/2, and 
finally ramped back to zero during the last quarter. 
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where v„i — 3d/{2tf) is the maximum trap velocity dur- 
ing the transport compatible with qoitf) = dm this 
scheme. Solving Eq. ^ for the previous qo{t) with initial 
conditions <7c(0) — (7c (0) = 0, and imposing continuity on 
qdt) and qdt) we find 

qc{tf)~qo{tf) = 9rf(l-2cos0)(sin2,/,)/(w2^2)^ 
9d 
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(smcf) + sm2(j) — sin3(^), (11) 
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FIG. 1: (Color online) (a) Displacement — qo versus time. 
Blue long dashed line: direct method, red short dashed 
line: inverse method (polynomial). Solid green line: inverse 
method + OCT. (b) Trap trajectories. Parameter values: 
d — 1.6 mm, tf = 20 ms, S ~ 0.162 mm, ujq — 2tt ■ 50 Hz. 



where (j) — a;oi//3. The final state of the transported 
BEC is given by Eqs. dZ]) and ([TT|). In general some 
excitation is produced, except for the discrete set of final 
times tf,N = 3(2iV+ l)7r/wo, iV = 0, 1,2,..., for which 
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and the transported state matches the eigenstate of the 
final Haniiltonian. The classically moving center of mass 
and the trap center stop at d, qdtf) = 0, qo{tf) = 0, with 
zero (classical) energy mqc{tf)^ /2 + mujQ[d — qdtf )]'^ /2 = 
0. Using this direct approach, the minimum final time 
which does not produce excitation is i/,o = Stt/wo {N = 
0). In our example, i/,o = 30 ms. For such short times 
the transport is not adiabatic. 

Thanks to the structure of the solution (O, we may 
apply a generalized inverse enginee ring method similar 
to the one for the linear case [1, 0, [Hf . The idea is to 
design qdt) first and deduce the transport protocol from 
it. We impose the conditions (|10|) and (fT2|) at i = and 
tf, and interpolate qc with a function, e.g. a polynomial 
with enough parameters to satisfy all these conditions. 
Then qo{t) is calculated via Eq. (|9]). An example is 
shown in Fig. 1 where we have chosen tf = 20 ms < 
tffi. By construction no final excitation is produced, and 
the final fidelity (overlap between the transported state 
and the ground state at tf) is one. Contrast this to the 
direct approach which, for tf — 20 ms, produces more 
transient excitation and a final excited state with nearly 
zero fidelity. 

In principle there is no lower limit to tf with the in- 
verse method, but in practice there are some limitations 



y. Smaller values oi tf increase the distance from the 
condensate to the trap center, see Eq. pTt . and the ef- 
fect of anharmonicity. There could be also geometrical 
constraints: for short tf, qo{t) could exceed the inter- 
val [0,d]. For the polynomial ansatz this happens [1] at 
tf — 2.505/a;o, t/ « 8 ms for the parameters of the ex- 
ample. Optical Control Theory (OCT) combined with 
the inverse method, see below, provides a way to design 
trajectories taking these restrictions into account. 

Anharmonic Transport. — The inverse method can also 
be applied to anharmonic transport by means of a com- 
pensating force 0]. To this aim, we consider a generic 
potential U{q — qo) and set a(t) = qo{t), uj{t) ^ luq = 0, 
/ = 0, and F — mqo in Eq. ([T]), so the CPE becomes 



.^ dijiq,t) 
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mqoq + U{q- qo) 



and the auxiliary equations ^ and ^ are satisfied triv- 
ially. Here we impose qo{0) = 0, qo{0) — 0, qo{tf) — 
d, qo{tf) = 0. We may optionally impose also ijo = 
at t = and tf. The function that must be interpolated 
is now qo(t), and again we consider a polynomial. For an 
arbitrary trap and tf = 20 ms, the maximal compensat- 
ing acceleration would be 23.1 m/s^. 

Optimal control theory. — Given the freedom left by the 
inverse method it is natural to combine it with OCT and 
design the trajectory according to relevant physical cri- 



teria 



21| . For harmonic transport, we have imposed the 
boundary conditions (fTO|) and (fT2| at t = and tf, but 
qo{t), and the polynomial ansatz for qc(t) are quite arbi- 
trary. As an example of the possibilities of OCT suppose 
that we wish to limit the deviation of the condensate from 
the trap center according to —S <qc — qol£^iS>(^ f^nd 
find the minimal time t f . The transport process given by 
Eqs. ([7]), (|T0| and (fT2|) can be rewritten as a minimum- 
time optimal control problem defining the state variables 
xi{t) and X2{t) and the control u{t), 

xi = a, X2 — a, u{t) = a — qo. (13) 

Equation (j9|) is transformed into a system of equations, 

X1—X2, ±2+^0^ = 0. (14) 

The OCT problem is to find -S < u{t) < S with 
u{0) = ultf) = 0, {a;i(0),a;2(0)} = {0,0}, and 
{xi{tf), X2{tf)} = {d,0} in the minimum final time tf. 
The optimal control Hamiltonian (23 | is He = Pia;2 — 
P2i^o^, where pi and p2 are conjugate variables. The 
Pontryagin maximality principle 22] tells us that for u{t), 
x{t) to be time-optimal, it is necessary that there exists a 
nonzero, continuous vector p(t) such that x = dHc/dp, 
p = —dHc/dx at any instant, the value of the control 
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FIG. 2: (Color online) Average fidelity of harmonic transport 
versus noise amplitude A. For tf — 20 ms: Qi/h = 0.05 m/s 
(solid, red line), Qi/fi ~ O.lm/s (dashed, blue line), Qi/h = 
0.2m/s (dotted, green line); for tf = 10 ms: Qi/h = 0.05m/s 
(red crosses), g^/h — O.lm/s (blue boxes), g^/h = 0.2m/s 
(green circles), uo = 2tt ■ 50 Hz. 



maximizes He, and Hc{p{t),x{t),u{t)) ~ c>0, with c a 
constant. The solution is of bang-bang type [23| . 



u{t) 



0, t = 

-6, <t<ti 

S, ti <t <tf 

0, t^tf 



where the initial and final discontinuities are chosen to 
satisfy the boundary conditions. Solving the system 
and imposing continuity on xi and X2 one finds for the 
switching and final times ti = t//2, tf = 2{d/Sy/^/ujo. 
The trap trajectory is deduced from Eq. ([1 



qoit) = 



0, t = o 

{l+uj^t^/2)6, 0<t<ti 
-[ujl{t-tf)^/2 + l]S + d, ti<t<tf ■ 

In Fig. [T] the displacement of the center of mass with 
respect to the trap center and the trap trajectory are 
plotted for this optimal trajectory. We have chosen 5 ~ 
0.162 mm so that the minimal final time is tf = 20 ms 
as in the previous example. 

Another important constraint might be that the center 
of the physical trap stays inside a given range (e.g. inside 
the vacuum chamber), i.e. the constraint is then < 
Qo{t) < Qt- Following the OCT procedure we finally get 



qo{t) 



0, t^o 

q^, <t <ti 

qi, ti <t <tf 

d, tf <t 



where wq^i 



1 



ujotf = uJoti + 



arccos ■ 



'{kld-^- qiiqi - qt)]/[id - qi)iqi - qt)]} ■ 
Noise. — In the following we investigate the effect of 
noise in harmonic transport. We assume that the cen- 
ter of the physical trap is randomly perturbed by the 
shift XC{t) with respect to qo{t). For the shifted trap 
center, Eq. ^ can be solved using the ansatz qc{t) = 
qc{t) + XPit) so that Pit) = J^"* drCir) sin (uJot-T), 



and ${t) = 1^0 /o " dr ({t) cos (wo< 
lution still given by Eq. (O. 
tp is independent of the chosen 



r), with the so- 
The fidelity at 



qc 



and d, F, 



C 



We assume 



/ dq exp (^ifXp{tp)q) x*{q + XI3{tF))x{q) 
now that (^(t) is white Gaussian noise, and average the 
fidelity over different realizations of C(t). The result 
can be seen in Fig. [2] for three values of g-^ and two final 
times, tf =20 ms and tf = 10 ms. The fidelity increases 
for smaller couplings g-^ and for the shorter time. 

Outlook — The above results may be extended to other 
physically motivated constraints, also to non-spherical 
traps with different frequencies uJx, Wy, lOz, rotations, 
and launching/stopping condensates up to/from a deter- 
mined velocity. 
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